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Abstract: Climate models predict that tropical lower-stratospheric humidity will increase as the 
climate warms. We examine this trend in two state-of-the-art chemistry-climate models. Under 
high greenhouse gas emissions scenarios, the stratospheric entry value of water vapor increases 
by ~1 part per million by volume (ppmv) over this century in both models. We show with 
trajectory runs driven by model meteorological fields that the warming tropical tropopause layer 
(TTL) explains 50-80% of this increase. The remainder is a consequence of trends in 
evaporation of ice convectively lofted into the TTL and lower stratosphere. Our results further 
show that, within the models we examined, ice lofting is primarily important on long time scales 
— on interannual time scales, TTL temperature variations explain most of the variations in lower 
stratospheric humidity. Assessing the ability of models to realistically represent ice-lofting 


processes should be a high priority in the modeling community. 
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Introduction 

Air traveling from the tropical troposphere into the tropical stratosphere transits the tropical 
tropopause layer (TTL) [Sherwood and Dessler, 2000], and the processes within this region 
provide primary control over the water vapor content of the stratosphere. In the following, we 
will refer to the water vapor mixing ratio of this air as HOentry. Over the past two decades, it has 
become generally accepted that H2Oentry variability is controlled by TTL temperature variability 
[e.g., Fueglistaler et al., 2009; Mote et al., 1996; Randel et al., 2004; Fueglistaler et al., 2005; 
Dessler et al., 2014; Wang et al., 2015]. This view posits that the TTL acts like a “cold trap,” 
where the humidity of lower stratospheric air is determined by the coldest temperatures 


experienced by the air as it crossed the TTL. 


Climate models have long predicted that H2Oecntry will increase over the next century [Gettelman 
et al., 2010; Kim et al., 2013], with important climatic [Forster and Shine, 1999; Solomon et al., 
2010; Maycock et al., 2013; Dessler et al., 2013] and chemical [Kirk-Davidoff et al., 1999] 
impacts. Despite the importance of these model results, few papers have analyzed the 
mechanism behind the overall increase in H2Oentry. Most papers that do view the problem 
qualitatively, finding that the increase in H2Oentry is roughly consistent with the long-term 
warming of the TTL [e.g., Fueglistaler and Haynes, 2005; Oman et al., 2008; Gettelman et al., 
2009; Garfinkel et al., 2013]. 


In this paper, we use a trajectory model driven by meteorology taken from climate models to 
quantitatively evaluate how much of the model trend in H2Oentry is due to changes in TTL 
temperatures and how much is due to water transport by other processes. We find strong 
evidence that while much of the future trend is due to a warming TTL, a significant fraction is 


due to increased transport of water in the form of convectively lofted ice. 


Models 

We analyze simulations from two chemistry-climate models (CCMs). These are similar to 
general circulation models, but with a more realistic stratosphere and higher vertical resolution in 
the TTL. As such, we expect CCMs to do a better job simulating H2Oentry than general 


circulation models. 
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GEOSCCM 

The Goddard Earth Observing System Chemistry Climate Model (GEOSCCM) couples the 
GEOS-5 general circulation model [Rienecker et al., 2008; Molod et al., 2012] toa 
comprehensive stratospheric chemistry module. The simulation used in this study has horizontal 
resolution of 2° latitude and 2.5° longitude with 72 vertical layers up to 0.01 hPa (80 km), with 
vertical resolution in the TTL of ~1 km. For our estimate of the GEOSCCM’s H2Oentry, we use 
the tropical average (30°N-30°S) 85-hPa volume mixing ratios. Averaging over 20°N-20°S 


yields nearly indistinguishable results. 


Prior versions of GEOSCCM have been extensively evaluated as part of the Chemistry-Climate 
Model Validation 1 (CCMVal-1) [Eyring et al., 2006] and CCMVal-2 [SPARC CCMVal, 2010], 
as well as in many other analyses [Strahan et al., 2011; Douglass et al., 2012; Oman and 
Douglass, 2014]. In this paper, we use a simulation from 1998-2099 driven by the RCP6.0 
scenario for greenhouse gases [van Vuuren et al., 2011] and the Al scenario for ozone depleting 
substances [World Meteorological Organization, 2011]. Sea surface temperatures and sea ice 
concentrations were prescribed from a CMIP5 simulation using the Community Earth System 


Model version 1 [Gent et al., 2011]. 


WACCM 

The Whole Atmosphere Community Climate Model (WACCM) is one of the available 
atmospheric components of the National Center for Atmospheric Research (NCAR) Community 
Earth System Model (CESM). WACCM includes processes essential to the simulation of the 
middle atmosphere such as nonlocal thermodynamic equilibrium radiative transfer, a non- 
orographic gravity wave drag parameterization, and a full representation of middle atmospheric 
chemistry that is coupled with radiation and dynamics [Hurrell et al., 2013; Marsh et al., 2013]. 
The simulation used here is a specified-chemistry version of WACCM (SC-WACCM) where the 
concentrations of radiatively/chemically active trace gasses are specified from existing WACCM 
simulations with interactive chemistry [Smith et al., 2014]. SC-WACCM was run at a horizontal 
resolution of 1.9° x 2.5° over 1955-2100 with the RCP 8.5 greenhouse gas scenario [van Vuuren 


et al., 2011]. This is a higher emissions scenario than that used in the GEOSCCM run, although 
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the effect of this on the analysis seems minor. The WACCM simulation includes a fully coupled 
ocean, land surface, and sea ice model as the other CESM components. For our estimate of the 
WACCM’s HoOentry, We use the same definition as for the GEOSCCM: tropical average (30°N- 


30°S) 85-hPa volume mixing ratios. 


The trajectory model 

We will compare estimates of H2Oentry from the CCMs to estimates from a domain-filling 
forward trajectory model [Schoeberl and Dessler, 2011]. In the version of the model analyzed 
here, an ensemble of 1350 parcels is initialized every day on an equal-area grid running from 
60°S to 60°N. The parcels are initialized at 370-K potential temperature (~16 km), which is 
above the level of zero net radiative heating in the tropics (~355-360 K) but below the tropical 
tropopause (~375-380 K). Each parcel is run forward until the parcel descends back into the 
troposphere, defined as pressures higher than 250 hPa (~10 km). All trajectory model runs 
include production of water vapor via methane oxidation, but that process is unimportant in the 


tropical lower stratosphere. 


The model uses the Bowman trajectory code [Bowman, 1993; Bowman and Carrie, 2002] to 
advect parcels, driven by 6-hourly instantaneous horizontal winds and 6-hourly average diabatic 
heating rates obtained from the GEOSCCM and WACCM runs. Each parcel is initialized with a 
water vapor mixing ratio of 200 parts per million by volume (ppmv). The mixing ratio is 
conserved along each trajectory, except when the relative humidity (RH) over ice of the parcel 
exceeds a pre-determined threshold [e.g., Schoeberl and Dessler, 2011], in this paper either 
100% or 80%. When parcels’ water vapor exceeds this threshold, the water vapor mixing ratio is 
instantly reduced until the RH equals the threshold value. The 100% threshold is frequently used 
in these types of analyses, but some CCMs begin dehydration below 100% [e.g., Molod, 2012], 
so this gives us some idea of the sensitivity of our results to differing thresholds. To estimate 
H2Oentry, we average the H2O mixing ratio of parcels between 79 and 93 hPa and between 30°N 
and 30°S. Dehydration events above 93 hPa do occur, but they remove relatively small amounts 


of water: the water vapor mixing ratio at 79 hPa is within a few percent of the value at 93 hPa. 
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We will refer to the model described in the previous paragraph as the 100% or 80% standard 
trajectory model, depending on the dehydration threshold. Despite the simplicity of this model, 
it has been shown to accurately reproduce many of the details of the water vapor distribution of 
the stratosphere [Schoeberl and Dessler, 2011]. Table 1 lists 21*-century average HzOentry in the 
standard trajectory models and the CCMs. The standard trajectory models do a good job of 
reproducing the CCMs’ value — to the extent they differ, the standard trajectory models tend to 
underestimate the CCMs. Most observational comparisons focus on water vapor anomalies 
(departures from the mean seasonal cycle), and the standard trajectory model does an excellent 
job reproducing observed anomalies [Schoeberl et al., 2012; Schoeberl et al., 2013; Dessler et al., 
2014; Wang et al., 2015]. 


CCM vs. trajectory model comparison 

The GEOSCCM predicts a change in H2Oenry over the 21" century (hereafter AH2O entry) of 0.87 
ppmvy, while the 100% and 80% standard trajectory model driven by GEOSCCM meteorology 
predicts AH2Oentry of 0.49 and 0.39 ppmv. The WACCM predicts AH2O entry of 1.09 ppmv, while 
the 100% and 80% standard trajectory models driven by WACCM meteorology predicts 
AH2Oentry of 0.86 and 0.70 ppmv. For all models, AH2Oentry is calculated as H2Oentry averaged 
over 2090-2100 minus the average over 2000-2010; values are also listed in Table 1. 


The disagreement between the CCMs and 100% standard trajectory model is shown graphically 
in Fig. 1. In the standard trajectory model, H2Oentry is entirely regulated by TTL temperature 
variations. The fact that the trajectory model mostly follows the GEOSCCM’s and WACCM’s 
H2Oentry variations lead us to our first main conclusion: TTL temperature variations are 
responsible for much of the trend in H»Oentry in the CCMs over the 21“ century. However, TTL 
temperature variations cannot explain all of the trends. In the GEOSCCM and WACCM, about 


50% and 20%, respectively, of the 21°-century trend must be due to other processes. 


A potential hint to explaining the discrepancy between the CCMs and the standard trajectory 
model is shown in Figure 2, which shows that convectively lofted ice-water content ((WC) in the 
GEOSCCM’s lower stratosphere increased significantly during the 21“ century. Convectively 
lofted IWC at 100 hPa more than doubles during the 21* century and increases by a factor of 
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about four at 85 hPa. The WACCM (not shown) only provides total IWC (the sum of convective 


and in situ ice) and that also shows an increase over the 21“ century. 


The convective injection of ice into the lower stratosphere, above the trajectories’ Lagrangian 
cold-point (LCP), where it can evaporate and moisten the stratosphere [e.g., Dessler et al., 2007; 
Schoeberl et al., 2014; Ueyama et al., 2015] may be the process missing from the standard 
trajectory model. LCPs in the 100% standard trajectory runs are found between 110 and 70 hPa, 


so the observations of convective ice at 100 and 85 hPa are consistent with this hypothesis. 


To test this idea, we run a second version of the trajectory model that includes the effects of 
convectively lofted ice, hereafter referred to as the “trajectorytice model”. In this model, we 
take the CCMs’ 6-hourly three-dimensional ice-water content (IWC) field and interpolate it onto 
each trajectory time step by linear interpolation in both time and space. At each time step, we 
assume complete evaporation of this ice into the parcel by adding the CCM’s IWC to the parcel’s 
water vapor, although we do not let parcels’ RH exceed the RH threshold, either 100% or 80%. 
Because we assume instant evaporation of the ice, we consider this to be an upper limit of the 


impact of convective ice evaporation on the water content of the TTL and lower stratosphere. 


Figure | shows that AH2Oentry from the 80% trajectory+ice model’s agrees more closely with the 
CCMs than either standard trajectory model (also seen in Table 1). The 100% trajectory+ice 
model (not shown) predicts slightly higher values of AH2Oecntry (Table 1). We noted above that 
the WACCM combines convective and in situ ice into one IWC variable, and we use that in the 
WACCM trajectorytice model. While this likely causes an overestimate of the evaporated ice in 
the WACCM.-based trajectory models, it may not be significant because in situ clouds tend to 
exist mainly in regions where RH is at or near saturation, so those clouds tend not to be 
evaporating. Table 1 also shows that the trajectory+ice models predict higher absolute values of 
H2Oentry than the CCMs, consistent with the idea that the trajectoryt+ice model is an upper limit 


on the effect of convective ice lofting. 


Figure 3 shows the spatial pattern of the change in HO mixing ratio at 100 hPa in the CCMs and 


two trajectory models over the 21“ century. It is clear that the trajectory+ice model more 
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accurately reproduces the spatial pattern found in both CCMs. The WACCM comparisons are of 
particular interest. For WACCM, the standard trajectory model actually does a reasonable job 
simulating the tropical average (e.g., Fig. 1 and Table 1), but Fig. 3 shows that it does a poor job 
simulating the spatial distribution of water. The trajectory+ice model, on the other hand, does a 
slightly better job simulating the tropical average, but a much better job reproducing the spatial 
distribution. The distribution at 85 hPa (not shown) also shows that the trajectoryt+ice model 


does a better job simulating the spatial distribution of HO. 


Are observations consistent with this result? 

We have demonstrated that convective ice lofting plays a key role in the long-term evolution of 
H2Oentry in the CCMs. One obvious question is whether observations are consistent with this. 
There have been many observational studies showing that convection penetrates into the tropical 
lower stratosphere [Alcala and Dessler, 2002; Dessler, 2002; Liu and Zipser, 2005; Dessler et al., 
2006; Rossow and Pearl, 2007], and there is also evidence that convective injection plays a role 
regulating the stratospheric water vapor budget [Moyer et al., 1996; Keith, 2000; Johnson et al., 
2001; Kuang et al., 2003; Hanisco et al., 2007; Corti et al., 2008; Khaykin et al., 2009; Schoeberl 
et al., 2014; Ueyama et al., 2015]. 


At the same time, many other analyses have concluded that observed H2Oentry variations over the 
last decade or so can be entirely explained by TTL temperature variations [e.g., Fueglistaler et 
al., 2009; Mote et al., 1996; Randel et al., 2004; Fueglistaler et al., 2005; Dessler et al., 2014; 
Wang et al., 2015]. This suggests a minor role for convective ice lofting, potentially 


contradicting results suggesting that convective lofting of ice is important. 


We can reconcile this seeming disparity by noting that observational studies necessarily cover 
short time periods. Over such short periods, the CCMs confirm that TTL temperature variations 
are indeed the main regulator of H2Oentry. This can be seen in Figure 4, which shows monthly 
H2O0entry anomalies from 2045-2055 from the CCMs agree with those from both the 100% 
standard trajectory model and the 80% trajectoryt+ice model. The clear message is that, while 
convective ice lofting is important for the long-term trend in H2Oentry in the CCMs, it does not 


play an important role in the CCMs’ short-term interannual variations. Thus, previous 
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conclusions that TTL temperature variability explains H2Oentry variability — based on a decade 
or so of data — should not be used to dismiss the potential importance of ice lofting in 21*- 


century trends. 


Nevertheless, the CCMs’ predictions of ice lofting into the lower stratosphere have not been 
quantitatively tested against observations. The CCMs’ predictions rely on their convective 
parameterizations, and until verified with observations, one can reasonably question the realism 
of their representation of the infrequent but intense convective systems that penetrate the 
stratosphere. In addition, the vertical resolution of the CCMs may not correctly resolve the top 
of convection, which could also bias the CCMs’ simulations. Validation of ice lofting in the 


CCMs should therefore be a high priority for the scientific community. 


Conclusions 
In this paper, we examine the long-term trend in H2Oentry, the humidity of air entering the tropical 
stratosphere, in two state-of-the-art chemistry-climate models (CCMs). The two models, the 


GEOSCCM and WACCM, both predict H2Oentry will increase over the 21“ century by ~1 ppmv. 


One hypothesis is that this trend is caused by a warming tropical tropopause layer (TTL). We test 
this by comparing H2Oentry from the CCM to that predicted by our trajectory models driven by 
the CCMs’ meteorology. The trajectory model sets water in each parcel to the minimum 
saturation mixing ratio the parcel experienced as it transited the TTL. We find that the warming 
of the TTL during the 21* century does indeed increase H2Oentry, but explains only 50-80% of the 
CCMs’ trends in H2Oentry. The remainder of the CCMs’ trends in H2Oentry must therefore be due 


to other processes. 


We identify the other process to be an increase in convectively lofted ice. If lofted above the 
Lagrangian cold point, the ice evaporates and moistens the stratosphere. Supporting this 
hypothesis is the fact that the CCMs predict increases in convectively lofted ice in the lower 
stratosphere. We tested the impact of this process by modifying the trajectory model to allow for 
the evaporation of convective ice. This trajectorytice model does a much better job simulating 


both the magnitude of the 21“ century trends and the spatial pattern. 
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We believe that solid evidence exists that trends in convectively lofted ice evaporation drives a 
significant part of the 21*-century trend in H2Oentry in the CCMs. This is mainly a long-term 
effect — on short time scales, the CCMs and trajectory models agree that TTL temperature 
variability drives most of the H2Oentry variability. This makes quantifying the impact of ice 
lofting in observational records difficult because observational records are generally too short for 
ice lofting to play a major role. Nevertheless, the importance of ice lofting on the long-term 
evolution of H2Oentry in CCMs should provide ample motivation to the community to study the 


fidelity of the CCMs’ representation of this process. 
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Figure captions 


Figure 1. Time series of AH2Oentry from (a) the GEOSCCM and two trajectory model runs driven 
by GEOSCCM meteorology and (b) from the WACCM and two trajectory model runs driven by 
WACCM meteorology. AH2Oentry is calculated by subtracting the average of the first 10 years 


from each time series. 


Figure 2. Annual and tropical average convectively lofted ice mixing ratio in parts per billion by 
volume (ppbv) from the GEOSCCM at 100 hPa (blue line, right-hand axis) and 85 hPa (red line, 
left-hand axis). 


Figure 3. The spatial distribution of the change in H2O over the 21 century at 100 hPa, 
calculated as the difference between the average of the last and first decades. Left column: 
GEOSCCM (top), GEOSCCM 80% trajectory+ice model (middle), GEOSCCM 100% standard 
trajectory model (bottom). Right column: the same quantities, but from WACCM. Each 


column’s color bar is located at the bottom of the column. 


Figure 4. Comparison between the CCMs, 100% standard trajectory model, and 80% 
trajectory+ice model over one decade (2045-2055). Quantities plotted are anomalies, which are 


the departures from that decade’s mean annual cycle. 


401 Table 1. Water vapor comparison between CCMs and trajectory model runs. The first column is 
402  H2Ocntry averaged over the 21" century. The second column is (AH2Oenty) is the change in 

403  H2Oentry over the 21“ century. The trajectory model listed under GEOSCCM use GEOSCCM 
404 meteorology while those listed under WACCM use WACCM meteorology. 


Model 21s'-century AH2Oentry 
avg. H2O0entry (ppmv) 
(ppmv) 
GEOSCCM 4.1 0.87 
100% standard trajectory 4.2 0.49 
80% standard trajectory 3.3 0.39 
100% trajectorytice 5.8 1.14 
80% trajectoryt+ice 4,7 0.92 
WACCM 4.7 1.09 
100% standard trajectory 4.0 0.86 
80% standard trajectory Be 0.70 
100% trajectorytice 6.5 1.20 
80% trajectoryt+ice 5.2 0.98 


405 
406 


407 
408 


409 
410 
411 


AH20 (ppmv) 


5 l 
2000 2020 2040 2060 2080 2100 


AH20 (ppmv) 


3 : b. WACCM 
- 2000 2020 2040 2060 2080 2100 


Figure 1. Time series of AH2O entry from (a) the GEOSCCM and two trajectory model runs driven 
by GEOSCCM meteorology and (b) from the WACCM and two trajectory model runs driven by 
WACCM meteorology. AH2Oentry is calculated by subtracting the average of the first 10 years 


from each time series. 
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Figure 2. Annual and tropical average convectively lofted ice mixing ratio in parts per 
billion by volume (ppbv) from the GEOSCCM at 100 hPa (blue line, right-hand axis) 
and 85 hPa (red line, left-hand axis). 
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Figure 3. The spatial distribution of the change in H2O over the 21“ century at 100 hPa, 
calculated as the difference between the average of the last and first decades. Left column: 
GEOSCCM (top), GEOSCCM 80% trajectory+ice model (middle), GEOSCCM 100% standard 
trajectory model (bottom). Right column: the same quantities, but from WACCM. Each 


column’s color bar is located at the bottom of the column. 
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425 
426 Figure 4. Comparison between the CCMs, 100% standard trajectory model, and 80% 


427  trajectory+ice model over one decade (2045-2055). Quantities plotted are anomalies, which are 
428 the departures from that decade’s mean annual cycle. 
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